The 1 < 2; < 5 Infrared Luminosity Function of Type I Quasars 
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ABSTRACT 

We determine the rest-frame 8/xm luminosity function of type I quasars over the redshift 
range 1 < z < 5. Our sample consists of 292 24/im sources brighter than 1 mJy selected from 
7.17 deg^ of the Spitzer Space Telescope MIPS survey of the NOAO Deep Wide-Field Survey 
Bootes field. The AGN and Galaxy Evolution Survey (AGES) has measured redshifts for 270 of 
the R < 21.7 sources and we estimate that the contamination of the remaining 22 sources by stars 
and galaxies is low. We are able to select quasars missed by ultra-violet excess quasar surveys, 
including reddened type I quasars and 2.2 < z < 3.0 quasars with optical colors similar to main 
sequence stars. We find reddened type I quasars comprise ~ 20% of the type I quasar population. 
Nonetheless, the shape, normalization, and evolution of the rest-frame 8/um luminosity function 
is comparable to that of quasars selected from optical surveys. The 8 fim luminosity function of 
type I quasars is well approximated by a power-law with index —2.75 ±0.14. We directly measure 
the peak of the quasar space density to be at z = 2.6 ± 0.3. 
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1. INTRODUCTION 

The luminosity function of Active Galactic Nu- 
clei (AGNs) is one of the principal constraints on 
the growth of supermassive black holes over the 
age of the Universe (e.g., Kauffmann & Haehnelt 
2000). Ideally, we would measure the bolomet- 
ric luminosity function of all AGNs as a func- 
tion of redshift. This is not trivial, as it requires 
deep multiwavelength data spanning X-ray to sub- 
mm wavelengths and sampling large comoving vol- 
umes. In practice, surveys select subsets of AGNs 
in relatively narrow wavelength ranges. For exam- 
ple, optical quasar surveys generally select quasars 
with blue spectral energy distributions and broad 
emission lines. While using a narrow wavelength 
range as a proxy for bolometric luminosity would 
appear risky, this approach seems to have been 
remarkably successful. Optical, radio, and X-ray 
luminosity functions of AGNs infer a z ~ 2.5 peak 
in the space density of the most luminous quasars 
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(e.g., Dunlop & Peacock 1990; Ueda at al. 2003; 
Croom et al. 2004). This agreement is reassur- 
ing, since the three selection methods have differ- 
ent biases - strong radio emission only occurs in 
a small subset of all AGNs, optical wavelengths 
are sensitive to dust absorption, and soft X-ray 
wavelengths are sensitive to gas absorption. 

It is the absorption of radiation by the gas and 
dust surrounding the black hole that creates the 
greatest systematic uncertainties. In particular, 
unified models of AGNs (e.g., Antonucci 1993) as- 
sume the presence of obscuring material with a 
non-isotropic spatial distribution that makes the 
observed spectrum a function of viewing angle. 
In unified models of AGNs, the obscuring mate- 
rial is often referred to as the torus, though the 
actual distribution of gas and dust may be sig- 
nificantly more complex and its properties could 
depend on luminosity, time or Eddington factors 
(e.g., Lawrence 1991; Hopkins et al. 2005). When 
AGNs are viewed from above the torus, optical 
continuum and broad emission lines associated 
with regions near the central black hole can be 
observed directly. When AGNs are viewed in 
the plane of the torus, optical continuum and the 
broad emission lines are absorbed. 

Selection biases due to absorption by dust or 
gas are minimized by selecting AGNs in the in- 
frared. In this paper we combine Spitzer Space 
Telescope 24 /im imaging of the Bootes field (Rieke 
et al. 2004; Houck et al. 2005) of the optical 
NO AO Deep Wide-Field Survey (NDWFS; Jan- 
nuzi & Dey 1999), with spectroscopy of quasar 
candidates by the AGN and Galaxy Evolution Sur- 
vey (AGES, Kochanek et al. in preparation) to 
determine the rest-frame 8 /im luminosity func- 
tion of type I quasars. This allows us to identify 
quasars that would be missed in optical surveys ei- 
ther because of extinction or because the quasars 
have optical colors typical of main sequence stars. 

The structure of this paper is as follows. In fj^] 
we describe the surveys used to produce the quasar 
sample. In Sj3]we characterize the sample, discuss 
its completeness and compare it to other selection 
methods available for the field. In ^ we derive 
the luminosity function of the type I quasars in 
the sample, and in we discuss the role of type 
II quasars. We summarize our results and discuss 
future prospects in |j6l Throughout this paper we 
define type I quasars as AGNs with optical and ul- 



traviolet emission lines broader than 1000 km 
and type II quasars as those with narrower lines. 
In this paper we adopt a cosmology of flm = 0.3, 
= 0.7, and Hq = 70 km s^^ Mpc. Throughout 
this paper we use AB magnitudes although the re- 
leased NDWFS, 2MASS, and USN0-A2 catalogs 
use Vega magnitude^ 

2. THE SURVEYS 

Our sample combines three surveys of Bootes - 
the NDWFS and Spitzer MIPS imaging surveys 
and the AGES spectroscopic survey. In this sec- 
tion we discuss the general properties of these 
three surveys and summarize the final multi- 
wavelength catalog. 

2.1. The NO AO Deep Wide-Field Survey 

The NDWFS is an optical and near-infrared 
imaging survey of two « 9.3 deg^ fields with 
NOAO telescopes (Jannuzi & Dey 1999). In this 
paper, we utilize the third data release of NDWFS 
BwRI imaging of the Bootes field, which is avail- 
able from the NOAO Science ArchivcH. A thor- 
ough description of the optical observing strategy 
and data reduction will be provided by Jannuzi et 
al. (in preparation). 

Object catalogs were generated using SExtrac- 
tor 2.3.2 (Bertin & Arnouts 1996) run in single- 
image mode. Detections in the different bands 
were matched if the centroids were within 1" of 
each other. For extended objects, detections in 
the different bands were matched if the centroid 
in one band was within an ellipse defined using 
the second order moments of the light distribu- 
tion of the object in another bancil. Throughout 
this paper we use SExtractor MAG_AUTO mag- 
nitudes (which are similar to Kron total magni- 
tudes; Kron 1980) for optical photometry, due to 
their small uncertainties and systematic errors at 
faint magnitudes. The 50% completeness limits 
of the SExtractor catalogs for point sources are 
Bw ^ 26.5, R ~ 25.5 and / ~ 25.3. 

We used SExtractor's star-galaxy classifier to 
classify objects as compact or extended. For 

^Bw,AB = Bw,Veg!L, RaB = Rvcga + 0.20, IaB = ^Vcga "t" 

0.45, and Jab = Jvegs. + 0.89 
^ http : / / www. archive .noao.edu/ ndwfs / 

^This ellipse was defined with the SExtractor parameters 
2 X A.WORLD, 2 x B.WORLD, and THETA.WORLD. 
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AGES target selection, objects were classified as 
compact if they had a CLASS_STAR parameter 
of > 0.8 in one or more optical bands. This re- 
sulted in a small number of emission line galaxies 
being erroneously classified as compact. For this 
paper, we classify objects as compact if they have 
a CLASS.STAR of > 0.7 in two or more of the 
three optical bands. Of the 288 24 fim sources 
with AGES spectroscopy which are classified as 
compact with our CLASS.STAR criterion, 20 are 
z < 1 galaxies which we reclassify as extended ob- 
jects for the remainder of the paper. Type I quasar 
candidates are selected from the sample of 24 /im 
sources with compact optical morphologies (f}3]). 
This reduces galaxy contamination while main- 
taining a high completeness for z > 1 quasars. 
Of the 1182 24 fim sources with R < 20.2 ex- 
tended counterparts and AGES spectroscopy, only 
two are z > 1 type I quasars. We do not attempt 
to measure the z < 1 type I quasar and Seyfert 
luminosity function, as the host galaxies of z < 1 
AGNs are sometimes detected and resolved in the 
deep NDWFS imaging resulting in a complex se- 
lection function that depends on host morphology, 
host luminosity, AGN optical luminosity, redshift 
and image quality. Our selection criteria also elim- 
inates type II quasars from the sample, which we 
will discuss in detail in a second paper. Our mea- 
surement of the 1 < z < 5 type I quasar luminosity 
function and our conclusions are not highly sensi- 
tive to our criteria for star-galaxy classification. 

2.2. The MIPS Survey 

The NDWFS Bootes field was observed simulta- 
neously at 24, 70, and 160 /im with the Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke 
et al. 2004) as part of the Spitzer IRS team's 
Guaranteed Time Observing programs (Houck 
et al. 2005). In this paper we only discuss the 
24 /xm data, as only five of the type I quasars in 
our sample were detected in the 70 /im and 160 /im 
imaging. These observations were performed with 
the "Medium Scan" technique. This observing 
mode is particularly suitable for high efficiency 
coverage of wide fields such as Bootes. The 24 /im 
detector uses a 2'.'55 per pixel array of 128x128 
elements with a delivered point source full-width 
at half maximum (FWHM) of ~ 6". The effective 
integration time per sky pixel was ~ 90 s. Data 
reduction was carried out using the MIPS Data 



Analysis Tool (Gordon et al. 2005). We corrected 
count rates for dark current, cosmic-rays, and flux 
non-linearities. These counts were then divided 
by flat -fields constructed for each position of the 
scan-mirror which is used in the "Medium Scan" 
observing mode. The images were corrected for 
geometric distortions and combined to produce 
the final 9.93 deg^ 24 /im mosaic. For this work, 
we selected 24 /im sources from the 8.22 deg^ of 
the 24 iim mosaic which overlaps the NDWFS and 
comprises of two or more individual scans. 

As almost all the sources in the MIPS mosaic 
are unresolved, source extraction and photome- 
try were performed using the stellar PSF fitting 
technique of the DAOPHOT software (Stetson 
1987). An empirical PSF was constructed from 
the brightest objects found in our mosaic, and it 
was subsequently fltted to all the sources detected 
in the map. Allowing for multiple-match fitting to 
deal with blended cases, we derived the flux den- 
sity of each source from the scaled fitted PSF and 
finally applied a slight correction to account for 
the finite size of the modeled point spread func- 
tion. Catalog depth was determined by adding 
artificial point sources to the MIPS mosaic and 
attempting to recover them with DAOPHOT. We 
estimate that the 80% completeness limit to be 
0.3 mJy. 

2.3. The AGN and Galaxy Evolution Sur- 
vey 

The AGN and Galaxy Evolution Survey (AGES, 
Kochanek et al. in preparation) is obtaining 
complete spectroscopic samples of galaxies and 
quasars selected over a broad range of wavelengths 
in the Bootes field of the NDWFS. The spectra 
are obtained with Hectospec, a 300-fiber robotic 
spectrograph on the MMT 6.5-m telescope (Fabri- 
cant et al. 1998; Roll, Fabricant, & McLeod 1998). 
Spectra were extracted and classified using two in- 
dependent reduction pipelines and then verified by 
visual inspection. One pipeline is based upon the 
SDSS spectroscopic data reduction software while 
the other is a set of customized IRAF scripts. The 
wavelength coverage of the spectra extends from 
3700A to 9200A, although some spectra are con- 
taminated redwards of 8500A by a red leak from 
an LED in the fiber positioner. 

The AGES spectroscopic targets were selected 
from the April 9, 2004 pre-release versions of the 
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NDWFS and MIPS 24 /xm source catalogs. How- 
ever, for the remainder of the paper we use the 
photometry of the third data release of the ND- 
WFS and July 5, 2005 pre-release version of the 
MIPS 24 fim catalog. These have less artifacts 
and slightly better photometric calibrations than 
the earlier catalogs. The spectra are mostly from 
AGES internal release 1.1 with some additional 
spectra of i?-band selected targets from AGES in- 
ternal release 2.0. AGES targeted 24 /xm sources 
with fluxes of 5*24 > 1 mJy (or [24] < 16.40 where 
[24] is the AB apparent magnitude) which had 
optically extended R < 20.2 counterparts or op- 
tically compact 17.2 < R < 21.7 counterparts. 
For AGES target selection, objects were classified 
as compact if their SExtractor CLASS_STAR pa- 
rameter was > 0.8 in one or more of the Bw, 
R, or /-bands. To reduce the number of bright 
(i? < 18) galactic stars in the spectroscopic sam- 
ple, compact sources were not targetted if they had 
2MASS counterparts bluer than J - [24] = -3.51. 
The roughly power-law spectra of quasars and the 
black body spectra of stars make the infrared color 
differences of these two populations very obvious 
(e.g.. Stern et al. 2005). These criteria led to 1571 
spectroscopic targets, and accurate spectroscopic 
redshifts were obtained for 1449 of them. 

2.4. The Multiviravelength Sample 

The final, multiwavelength catalog consists of 
24 /im sources brighter than 1 mJy with good 
24 /Ltm and optical imaging inside the primary 
AGES target regions. For saturated objects in 
the NDWFS images, we have used the fourth data 
release of the Sloan Digital Sky Survey (SDSS 
DR4; Adelman-McCarthy 2005) for the photome- 
try and classifications. We excluded sources with 
R > 17.2 counterparts within (20.2 - R)" of 
R < 15.2 USN0-A2 stars (Monet et al. 1998) 
where our ability to detect faint optical counter- 
parts is compromised. We also excised > 1' radius 
regions around Hipparcos stars (Esa 1997), which 
are occasionally broken up into multiple spurious 
USN0-A2 catalog entries. This left us with an ef- 
fective survey area of 7.17 deg^. Objects with spu- 
rious R > 21.7 counterparts (e.g., edges of i? 18 
star halos) have been removed from the final cata- 
log. The sky distribution of objects in the catalog 
is shown in Figure [TJ 

Table [1] presents the number counts of the 




Fig. 1. — The sky distribution of 24 /im sources 
brighter than 1 mJy in the multiwavelength cat- 
alog. Compact objects with and without spectra 
are shown with solid and open stars respectively. 
Gray circles denote the locations of sources with 
extended optical morphologies. The semicircular 
edges of the sample are an artifact of the circular 
1 degree field of view of the Hectospec instrument. 



24 nm sources and Figure [2] shows the distribution 
of their i?-band magnitudes. Most of the 24 /im 
sources are brighter than R = 22, and most of 
these sources are either quasars or z < 1 galaxies. 
The 17% of 24 /xm sources with R > 21.7 coun- 
terparts include include ultra-luminous infrared 
galaxies and obscured quasars at z > 1 (Houck 
etal. 2005; Martfnez-Sansigreet al. 2005). Type I 
quasars have relatively bright optical counterparts 
and mostly compact optical morphologies. AGES 
spectroscopic redshifts are available for 63% of the 
524 > 1 mJy sources and 92% of the 24 /xm sources 
satisfying the AGES selection criteria. 

3. THE TYPE I QUASAR SAMPLE 

We define our type I quasar sample to be 24 /xm 
sources brighter than 1 mJy with 17.2 < R < 23.5 
compact optical counterparts. The compact op- 
tical morphology is due to the observed optical 
emission being dominated by a relatively unob- 
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Fig. 2. — The distribution of i?-band magnitudes for all {left) and optically compact (right) 24 /j,m. sources. 
The distribution of objects with spectroscopic redshifts (from both the literature and AGES) is shown with 
the shaded histogram. Most (86%) of objects with compact optical morphologies have AGES spectroscopic 
redshifts. 



scured AGN. Our sample is incomplete for z < 1 
quasars and Seyferts, where AGN host galaxies are 
sometimes resolved in the NDWFS imaging and 
AGN host galaxies can contribute significantly to 
the observed 24 /xm flux. The bluest spectroscopic 
quasar in our sample has a color of i? — [24] = 2.5 
(where R — [24] is the apparent color in AB mag- 
nitudes), so we have excluded R— [24] < 2 sources 
from the sample with the expectation that this will 
not effect the sample completeness for highly lumi- 
nous quasars. Of the 47 objects excluded with this 
color criterion, 90% are brighter than R = 17 and 
are almost certainly stars. AGES spectroscopy 
identified 20 z < 1 galaxies which were erroneously 
classified as compact in the NDWFS, and they 
were reclassifed as extended and excluded from the 
quasar sample. Two R < 20.2 2 > 1 quasars which 
were classified as extended in the optical imaging 
have been added to the quasar sample. The fi- 
nal sample contains 292 quasars and quasar can- 
didates. Figure[3]shows the redshift distribution of 
the 270 objects in the sample with AGES spectro- 
scopic redshifts. All the quasars without redshifts 
are at z < 3.8 (or are blends with z < 3.8 objects), 
as they all have Bw detections indicating that the 



Lyman limit is within or blueward of our Byy fil- 
ter. The vast majority (83%) of the quasars with 
spectra and 59% of the candidates without spectra 
have one or more photons in the XBootes Chandra 
imaging survey (Murray et al. 2005). In contrast 
only 23% of the z < 1 galaxies with S24 > 1 mJy 
are detected in XBootes. Quasar candidates with- 
out spectra have the same likefihood (~ 10%) of 
having 1.4 GHz counterparts (from de Vries et al. 
2002) as spectroscopically confirmed quasars. We 
therefore expect most of the 22 quasar candidates 
without spectra to be quasars. 

While our i?-band spectroscopic flux limit is 
less sensitive to extinction than standard ultravi- 
olet excess (UVX) searches, it still corresponds to 
the rest-frame ultraviolet emission of the quasars. 
As shown in Figure IH the apparent colors and 
magnitudes of quasars show some evidence for ob- 
scuration since fainter quasars tend to have redder 
colors that could be explained by modest amounts 
of extinction. This is consistent with modest dust 
reddening producing the broad red tail of type I 
quasar colors (e.g., Richards et al. 2001). 
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Spectroscopic Redshift 



Fig. 3. — The spectroscopic redshift histogram of type I quasars. The sample is incomplete at z < 1 as 
some quasars are classified as extended in the deep NDWFS imaging. The curve is derived from the best-fit 
luminosity function, and is a good approximation to the z > 1 data. 




Fig. 4. — The apparent colors of type I quasars. The symbol size and color are a function of redshift. 
Objects with spectroscopic and photometric redshifts are shown with solid and open symbols respectively. 
Arrows denote the apparent color vectors for z = 2 dust extinction with rest-frame E{B — V) — 0.1, using 
the approximation of Calzetti (1999). Only objects with unsaturated Bw photometry are shown in the left 
panel. Relatively little dust obscuration can explain much of the dispersion of type I quasar colors. 
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Fig. 5. — SDSS DR4 color-magnitude (left) and color-color (right) diagrams for the 24 /im selected type I 
quasar sample. As in Figure IH symbol size and color are a function of redshift, and open symbols denote 
objects without spectroscopic redshifts. In the right panel, the 2SLAQ selection criteria for g < 21.15 and 
g < 21.85 are shown with bold and narrow lines respectively. Quasars selected by the 2SLAQ criteria are 
shown with squares. Of the 142 1.0 < z < 2.2 quasars in the 24 /im sample, 112 (79%) are also selected by 
the UVX criteria of 2SLAQ. 



3.1. Comparison with UVX Surveys 

As a check of our sample completeness and to 
compare our selection criteria with traditional op- 
tical surveys, we compared our sample to a sam- 
ple of UVX quasar candidates. These were se- 
lected using the criteria of the 2dF-SDSS Lumi- 
nous Red Galaxy and Quasar Survey (2SLAQ; 
Richards et al. 2005) and SDSS DR4 imaging cat- 
alogs, which overlap almost all of the Bootes field. 
The 2SLAQ criteria do not include morphology 
cuts for most objects, but do include optical color 
cuts which limit the ability of 2SLAQ to detected 
reddened and z > 2.2 quasars. We restricted the 
comparison to i? < 19 quasars because at fainter 
magnitudes the blue quasars will tend to have 
S24 < 1 mJy (see Figure E]). Of the 46 i? < 19 
UVX candidates in Bootes, 38 (83%) are in our 
24//m quasar sample. The remaining i? < 19 
UVX candidates are starburst galaxies or z < 0.3 
quasars whose host galaxies are detected in the 
NDWFS imaging (and are classified as extended 
sources). Our 24/im quasar sample is therefore 
highly complete for type I quasars selected with 
the UVX technique, particularly at z > 1 where 



quasar hosts are unresolved in the NDWFS imag- 
ing. 

We now consider the completeness of UVX sur- 
veys for quasars selected via their 24 /im flux. Of 
the 23 i? < 19 quasars with 1.0 < 2: < 2.2 in 
our 24 /im sample, 19 (83%) are selected by the 
2SLAQ UVX criteria. The four 1.0 < z < 2.2 
quasars not selected by the 2SLAQ criteria all have 
red colors {u- g > 0.5 and R - [24] > 3.7). Given 
the problems of small number statistics, this is 
consistent with the estimate of 90% completeness 
for the UVX selected quasars from the 2QZ cat- 
alog by Smith et al. (2005). However, as shown 
in Figure [SJ the completeness of the UVX sam- 
ples of quasars above a fixed infrared flux is lower. 
This is relevant for understanding the accretion 
history of quasars as the infrared is almost cer- 
tainly more strongly correlated with the bolomet- 
ric luminosity than the ultraviolet. If we select 
the 142 1.0 < z < 2.2 quasars brighter than 
1 mJy at 24 fim, we find 79% are selected by 
the 2SLAQ criteria. The quasars not selected by 
the 2SLAQ criteria are sometimes optically blue 
and faint, but generally they are red. The me- 
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dian apparent colors of the 30 1.0 < z < 2.2 not 
selected by the 2SLAQ criteria are u — g = 0.66 
and R — [24] — 4.62. Clearly, there is a signif- 
icant population of bolometrically luminous red 
type I quasars which are mostly undetected by tra- 
ditional UVX quasar surveys. Red type I quasars 
comprise ^ 20% of the type I quasar population. 

3.2. Morphological Completeness 

We can also check the consequences of includ- 
ing only compact sources in our sample. This cri- 
terion certainly removes type II quasars from the 
sample, as their optical emission is dominated by 
their host galaxies, but it will also remove z < 1 
type I quasars from our sample when their host 
galaxies are detected in the deep NDWFS images. 
Thus, when we evaluate the quasar luminosity 
function in SJH we restrict the analysis to z > 1 
quasars. AGES targeted both compact and ex- 
tended 17.2 < R < 20.2 sources, so for this magni- 
tude range we can estimate the incompleteness due 
to compact selection criteria. Of 130 24/im sources 
with z > 1 and R < 20.2, only two quasars (1.5%) 
were classified as extended sources. These objects 
have been added to our spectroscopic quasar sam- 
ple. One is close to a foreground star while the 
other is a blend with a galaxy, which may or may 
not be physically associated with the quasar. At 
20.2 < R < 21.7, we assume the z > 1 quasar sam- 
ple is 1.5% incomplete due to quasars being clas- 
sified as extended in the optical imaging. Thus, 
the morphology cut should only affect the sample 
completeness at levels well below our statistical 
uncertainties. 

3.3. Spectroscopic Completeness 

We can see in Figures [2] and [4] that the vast 
majority of S'24 > 1 mJy type I quasars are 
brighter than our spectroscopic limit of R = 21.7 
and have spectroscopic redshifts. The number of 
S'24 > 1 rnJy type I quasars is clearly decreasing 
with increasing magnitude at i? > 21.7. Unlike 
UVX criteria, our selection criteria do not intrude 
into the main locus of type I quasar colors. We can 
use a relatively simple model of the spectroscopic 
completeness, and our estimate of the luminosity 
function is not particularly sensitive to changes in 
this model. 

To model the spectroscopic completeness, we 



assume the distribution of rest-frame R — [24] col- 
ors for z > 1 type I quasars is neither a function of 
luminosity or redshift. We measure this distribu- 
tion using the observed colors of z > 1 quasars 
and /c-corrections derived from the higher mid- 
infrared (MIR) flux quasar spectral energy dis- 
tribution (SED) of Hatziminaoglou et al. (2005). 
When spectroscopic redshifts are unavailable, we 
estimate photometric redshifts by searching for 
the nearest neighbor in color-space with a spectro- 
scopic redshift. As shown in Figure [6l the appar- 
ent R— [24] colors of quasars are a weak function of 
redshift so relatively inaccurate photometric red- 
shifts can be used to determine rest-frame colors 
(but not luminosities). 

Type I quasars fainter than R = 23.5 are not in- 
cluded in our sample, as they are not classified as 
compact sources. We include a completeness cor- 
rection for these objects by assuming their color 
distrbution is identical to quasars with higher 
24 /ini fluxes and that the shape of their 24 /im 
number count distribution is identical to that of 
bluer quasars. We estimate there are four z > 1 
type I quasars fainter than R = 23.5, and include 
this in our completeness estimates though this has 
little effect on our measured luminosity function. 

Our estimate of the rest-frame R— [24] distribu- 
tion for z > 1 quasars is shown in the left panel of 
Figurc[7l To derive the distribution of apparent R- 
band magnitudes as a function of both 24 /im flux 
and redshift, we used this model and applied the 
fc-corrections derived from the higher MIR quasar 
SED of Hatziminaoglou et al. (2005). We mul- 
tiplied the distribution of i?-band magnitudes by 
the measured completeness of as a function of R- 
band magnitude to determine the spectroscopic 
completeness of the sample as a function of 24 /im 
flux and redshift. As shown in the right panel of 
Figure [7l the spectroscopic completeness peaks at 
a 24 fiia flux of 7m Jy, which corresponds to 
17.2 < R < 20.2 where the AGES spectroscopic 
completeness plateaus at ~ 95%. The spectro- 
scopic completeness then declines to ~ 80% at 
both 1 mJy (because of the unmeasured redshifts) 
and at 15 mJy (because of the AGES R > 17.2 
criterion). 
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Spectroscopic Redshift Spectroscopic Redshlft 

Fig. 6. — The apparent colors of 24 fim sources in the Bootes field as a function of spectroscopic redshift. 
Optically compact objects are shown with stars while extended objects are shown with dots. Two quasar 
templates from Hatziminaoglou et al. (2005) are also shown. The dashed line denotes the template derived 
from 35 SDSS quasars in the SWIRE ELAIS-Nl field while the solid line denotes the higher MIR flux quasar 
template, which is a better fit to the median color of our sample. The apparent R — [24] color of type I 
quasars is a weak function of redshift in the range 1 < z < 5. The spectroscopic sample contains very red 
and bolometrically luminous quasars, though at the limits of the 24 fim and AGES spectroscopic samples 
we are restricted to i? — [24] < 5.3 objects. 




Rest-frame R-[24] (AB) 24nm flux (mJy) 



Fig. 7. — The distribution of rest-frame colors (left panel) and the estimated spectroscopic completeness 
(right panel) for the 24 /im selected type I quasar sample. For the redshift range 1 < z < 4, all S'24 > 1 mJy 
type I quasars with restframe R— [24] < 6 are brighter than our R — 21.7 spectroscopic limit. The plateau of 
the spectroscopic completeness corresponds to 17.2 < R < 20.2, where the AGES spectroscopic completeness 
for optically compact sources is ~ 95%. 
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4. LUMINOSITY FUNCTION OF TYPE 
I QUASARS 

The luminosity function of quasars is one of the 
principal constraints on the accretion history of 
the most massive blackholes (e.g., Kauffmann & 
Haehnelt 2000). At the highest luminosities the 
fraction of type II quasars is thought to be low 
(e.g., Ueda et al. 2003; Hao et al. 2005), although 
current studies do not include radio-quiet quasars 
obscured in both the optical and X-ray. If there 
are few type II quasars at very high luminosities 
then the type I quasar luminosity function pro- 
vides a good approximation of the overall quasar 
space density at the highest luminosities. 

We evaluate the luminosity function using 183 
z > 1 quasars with spectroscopic redshifts. As 
discussed earlier, the deep NDWFS images can 
sometimes detect the host galaxy of z < 1 quasars 
so modeling the completeness would require a de- 
tailed treatment of star/galaxy classification as 
a function of quasar luminosity, host luminosity, 
host morphology and image quality. Several z < 1 
extended objects are classified as broadline AGNs 
by the AGES pipeline and are not included in our 
sample of 292 type I quasars. The present sam- 
ple does not include z > 5 quasars because of the 
bright i?-band magnitude limit. In the 2005 AGES 
sample, which uses a fainter /-band selection cri- 
terion, appreciable numbers of z > 5 quasars are 
found, but they are fainter than S24 = 1 mJy. 

We estimated the rest-frame 8 fj,m absolute 
magnitudes of the sample using the AGES spec- 
troscopic redshifts, the 24 /im flux, and the higher 
MIR flux quasar SED of Hatziminaoglou et al. 
(2005). For the redshift range 1 < 2: < 5, ob- 
served 24 /j,m corresponds to rest-frame 12 /im to 
4 urn, so our 8 fiia absolute magnitudes are not 
very sensitive to our choice of quasar SED. Us- 
ing a bluer quasar SED template, derived from 
35 SDSS quasars in the SWIRE ELAIS-Nl field 
(Hatziminaoglou et al. 2005), alters the absolute 
magnitudes of individual quasars by only a few 
tenths of a magnitude, our luminosity function 
parameters by < la, and our conclusions little. 
The distribution of absolute magnitudes as a func- 
tion of redshift is shown in Figure [H Most objects 
have 8 iiui absolute magnitudes brighter than —25, 
which corresponds to 1.6 x lO'^^ergs s~^ in lyF^. 

We estimate the luminosity function using the 




Spectroscopic Redshift 



Fig. 8. — The distribution of the 8 /im absolute 
magnitudes of the type I quasar sample as a func- 
tion of spectroscopic redshift. Absolute magni- 
tudes were determined using the measured 24 /im 
fluxes, spectroscopic redshifts, and the higher MIR 
flux template of Hatziminaoglou et al. (2005). 
Since rest-frame 8 fim roughly corresponds to ob- 
served 24 /im, we have measured 8 /im absolute 
magnitudes and the 8 /im luminosity function. 
The Mb < —22.3 criterion for quasars (Veron- 
Cetty & Veron 2003) roughly corresponds to an 
8 /im AB absolute magnitude limit of —25.3. 

maximum likelihood method of Marshall et al. 
(1983). We assume the luminosity function is a 
power-law of the form 

ML,z)dL^ (z^) 

where a is a constant and L*{z) describes the 
pure luminosity evolution of the quasar luminosity 
function. This can also be written as a function of 
absolute magnitude; 

MM, z)dM cx ioO-4("+i)[M*(-)-M]^^j_ (2) 

Optical surveys generally require models where 
the slope of the luminosity function becomes shal- 
lower for low luminosity quasars (e.g., Schmidt 
& Green 1983; Boyle, Shanks, & Peterson 1988; 
Groom et al. 2004; Wolf et al. 2003b; Richards 
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Fig. 9. — The 8 /zm luminosity function of z = 2 type I quasars. The sohd hne is the best-fit power-law while 
the datapoints are determined using 1/Vmax with quasar luminosities adjusted for evolution. For clarity, we 
only show error bars and 3ct upper limits for the 1/Vmax luminosity function determined with 1 < z < 5 
quasars. For comparison, the dashed line shows the 2QZ luminosity function (Croom et al. 2004) where we 
have assumed that Bj — [8] = 3.08 (AB) for blue quasars. An AB absolute magnitude of —28.5 corresponds 
to 4.11 X 10'*^ergs s~^ in i^Fi,. 



et al. 2005), but a single power-law fits our data 
reasonably well. When we tried fitting a double 
power-law to the data, we were unable to precisely 
constrain the faint-end slope due to the limited 
size and luminosity range of our sample. We mod- 
eled the evolution of the quasar luminosity func- 
tion using pure luminosity evolution where 

L*{z) = L* {0)10''^'+''^'' (3) 

with fci, k2, and being constants. This is similar 
to the polynomial used by optical quasar surveys 
such as the 2QZ, and aids comparisons with the 
previous surveys. Unlike optical quasar samples, 
we have no difficulty finding quasars at 2.3 < z < 
3.0 (where the optical colors match those of main 
sequence stars), so we can track the evolution of 
the luminosity function through its peak. This 
forces us to add the k^z^ term to the polynomial 



so that the evolution on the high redshift side of 
the peak is not forced to be identical to that on 
the low redshift side of the peak by the symmetry 
of the equation. This pure luminosity evolution 
model is sufficient for a sample of 10^ quasars 
- larger samples containing ^ 10* quasars require 
more complex models (e.g., Croom et al. 2004). 

In Figure[n]we plot our best-fit luminosity func- 
tion for z = 2 type I quasars and our best-fit pa- 
rameters are summarized in Tablc[2] Only the 183 
quasars with spectroscopic redshifts of z > 1 were 
used to measure the luminosity function. We used 
the completeness estimates described in §3.21 and 
§3.31 to correct for the spectroscopic and morpho- 
logical incompleteness of the sample. We also gen- 
erated an estimate of the z = 2 luminosity func- 
tion using the 1/Vmax method after attempting to 
remove the effects of evolution. For each quasar. 
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(2005) template derived from 35 SDSS quasars approximates the median color for these samples. This is 
only an approximation since it ignores the spectral diversity of blue quasars. The space density of bright 
quasars rapidly evolves and peaks at2: = 2.5±0.3. 



we estimated a z = 2 absolute magnitude by as- 
suming 

M(2) ^ M{z) + 2.5{kiz + k2Z^ + k3Z^) 

-2.5{2ki + 4k2 + 8k3) (4) 

with our best- values of fci , k2, and fcs (Table [2]) . 
This allows us to use quasars from a broad range 
of redshifts to estimate a 1/Vmax luminosity func- 
tion. We have used a statistical representation of 
the evolution rather than a model for the evolu- 
tion of individual quasars. A large number of black 
holes undergoing relatively brief episodes of quasar 
activity can produce a luminosity function that 
mimics a smaller population of quasars undergoing 
luminosity evolution (e.g., Kauffmann & Haehnelt 
2000). Our measurements of the z — 2 quasar lu- 
minosity function using the l/Vmax method are 
plotted in Figure [5] and tabulated in Table [31 Our 
best fit power-law and l/Vmax luminosity func- 
tions agree well. The very faintest quasars show a 
slight departure from the power-law, but we can- 
not detect it with any precision when we include 



a break in our maximum likelihood models. Fig- 
ure O shows that the observed and model distribu- 
tion of the 1 < 2 < 5 quasars in redshift agree well, 
though our results cannot be naively extrapolated 
to other redshifts. 

Our best-fit parameters for the quasar luminos- 
ity function are summarized in Table [2] along with 
the luminosity function parameters for the 2QZ 
(Croom et al. 2004) for comparison. Our lumi- 
nosity function is not as steep as the bright-end 
of the 2QZ luminosity function. However, we are 
measuring quasars that are slightly brighter than 
the break in the 2QZ luminosity function, where 
the power-law index changes from —3.3 to ~ —1.5. 
As shown in Figure [H the z = 2 luminosity func- 
tions of this work and the 2QZ are similar for 
quasars within our luminosity range. Figure [10] 
plots the measured evolution of the quasar space 
density from our 24 /im quasar sample and the op- 
tically selected 2QZ and SDSS samples. The three 
surveys indicate a peak in quasar activity occurs 
at z ~ 2.5. 
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If the dust obscuration of AGNs is a function 
of bolometric luminosity (e.g., Hao et al. 2005), 
then the shape of the quasar luminosity function 
should vary with color. For example, if obscura- 
tion increases with decreasing luminosity, then the 
luminosity function of reddened quasars should be 
steeper than that of blue quasars. We divided our 
sample into quasars with R — [24] colors that are 
redder and bluer than the higher MIR flux quasar 
template of Hatziminaoglou et al. (2005). This 
rest-frame color cut oi R— [24] = 4.4 corresponds 
to an apparent color cut oi R — [24] ~ 3.8 over 
the redshift range 1 < z < 4. As the spectro- 
scopic subsamples contain ~ 90 quasars each, the 
uncertainties for a are prohibitively large if we let 
evolution parameters {ki, k2, and k^) float. We 
have therefore fixed the values of the evolution pa- 
rameters to the fit derived with the entire quasar 
sample. 

Our best fit luminosity function parameters are 
given in Tabled The power-law indices (a) of the 
blue and red subsamples differ by only la. This 
is not completely unexpected as the range of dust 
obscuration probed by our samples is limited. As 
shown in Figured at z = 2 our 17.2 < i? < 21.7 
spectroscopic sample spans a range of E{B — V) 
values of ^ 0.4. In contrast, the range of E{B — V) 
values spanned by Seyfert Is and lis is ^ 100 (e.g., 
Clavel et al. 2000). At present, the dependence 
of the dust obscuration on bolometric luminosity 
has only been seen in comparisons of type I and 
II AGNs (Hao et al. 2005), so a complete sam- 
ple of type II quasars may be required to under- 
stand the relationship between dust obscuration 
and luminosity for the most bolometrically lumi- 
nous quasars. 

5. TYPE II QUASARS 

Ideally we would measure the luminosity func- 
tions of both type I and II quasars. However, 
spectroscopic redshifts for z > 1 type II quasars 
are non-trivial as their optical counterparts can be 
extremely faint (e.g., Houck et al. 2005; Martfnez- 
Sansigre et al. 2005). Photometric redshifts are 
also difficult as the broad-band SED includes un- 
known contributions from the AGN, host galaxy, 
and dust (Zakamska et al. 2005). In this section, 
we provide a brief discussion of the type II quasar 
candidates. A spectroscopic campaign to measure 



the redshifts of the most luminous type II quasar 
candidates is ongoing, and a future paper will dis- 
cuss the luminosity function of these objects in 
detail. 

We designate all 24 fim sources brighter than 

I mjy with extended R > 21.7 optical counter- 
parts (or counterparts too faint for reliable mor- 
phological classification) as z > 1 type II quasar 
candidates. Atz^l,i? = 21.7 corresponds to 
a blue ~ L* galaxy (Wolf et al. 2003a). Unlike 
radio and X-ray surveys, we can detect z > 1 type 

II quasars which are both radio-quiet and Gomp- 
ton thick. The main contaminants of the sample 
will be starburst galaxies, z < 1 type II AGNs, and 
reddened type I quasars. There are 376 z > I type 
II quasar candidates in our survey area. Only 20 of 
these objects have spectroscopic redshifts, mostly 
from Houck et al. (2005), and all but one of these 
redshifts are in the range 0.7 < z < 3.0. 

The optical colors of the type II quasars candi- 
dates are shown in Figure [TlJ The locus of blue 
Bw — R and red R — I colors are consistent with 
1 < z < 3 host galaxies. Type II candidates 
with counterparts in the radio catalog of de Vries 
et al. (2002) and the XBootes Ghandra imaging 
survey (Murray et al. 2005) are also shown in 
Figure [11] As the 0.5 to 7 keV background of the 
XBootes survey is only ~ 0.03 counts within a 2" 
radius aperture, the detection of a single X-ray 
photon within this aperture is significant in the 
vast majority of cases, and at z = 1 a single photon 
in XBootes corresponds to an X-ray luminosity of 
~ lO''^ ergs s~^. For a z > 1 source this is a strong 
indicator of AGN activity, as the implied luminos- 
ity is an order of magnitude more than low red- 
shift starbursts including Arp 220, which has an 
observed X-ray luminosity of 1.4 x 10"'^ ergs s~^ 
(McDoweU et al. 2003). One-third of the type II 
quasar candidates have one or more X-ray photons 
in the Ghandra images, and this increases to 45% 
for candidates brighter than 2 mJy. 

We have determined upper limits for the frac- 
tion of z > 1 obscured quasars as a function of 
flux by dividing the number of z > 1 type II can- 
didates by the total number of z > 1 type I and 
type II candidates. Gontamination of our z > 1 
type I sample is low and contamination of the 
z > 1 type II candidate sample by z < 1 AGNs 
and starbursts could be high, so this fraction is an 
upper limit. Upper limits for the fraction of ob- 
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Fig. 11.— The BwRI colors of i? > 21.7 type 
II quasar candidates. Symbol size is a function 
of 24 /xm flux and objects with radio (bold sym- 
bols) and X-ray counterparts (crosses) are high- 
lighted. The optical colors of most objects are 
consistent with 1 < z < 3 host galaxies. Can- 
didates with X-ray counterparts are almost cer- 
tainly AGNs, and they comprise 32% and 45% of 
the candidates brighter than 1 mJy and 2 mJy 
respectively. 

scured quasars as a function of 24 fim flux are plot- 
ted Figure [121 This should not be confused with 
the fraction of obscured quasars as a function of 
bolometric luminosity, which would require accu- 
rate redshift information. We do not know if the 
obscured quasars have systematically lower red- 
shifts (and bolometric luminosities) than the type 
I quasars with comparable 24 /im fluxes. We also 
do not know the extent of dust absorption of rest- 
frame ^ 8 /im emission in these AGNs. Despite 
these caveats, it is reasonable to conclude that 
the ratio of obscured and unobscured AGNs near 
the break in optical luminosity function (roughly 
— 1 mJy) is of order unity. This is consis- 
tent with studies of X-ray selected AGNs, opti- 
cally selected Seyferts, and small samples of in- 
frared selected z > 1 type II quasars (e.g., Stef- 
fen et al. 2003; Ueda et al. 2003; Hao et al. 
2005; Martfnez-Sansigre et al. 2005). The frac- 
tion of obscured quasars decreases with increasing 



Fig. 12. — Upper limits on the fraction of ob- 
scured quasars as a function of observed 24/im 
flux. The shaded area denotes the region below 
our best estimate of the obscured fraction plus the 
1(7 statistical uncertainties. Our sample of type II 
quasar candidates could be contaminated by star- 
burst galaxies, so we can only set upper limits in 
the absence of spectroscopic confirmation. The 
fraction of obscured quasars could be as high as 
70% at 1 mJy, and then decreases with increasing 
24 iim flux. 

flux, which is similar to the trend observed as a 
function of luminosity in X-ray and optical sam- 
ples (e.g., Ueda et al. 2003; Hao et al. 2005). If 
the obscured fraction is decreasing rapidly with in- 
creasing luminosity, the type I quasar luminosity 
function approximates the luminosity function of 
all quasars at the highest luminosities. However, 
to confirm this will require complete spectroscopic 
samples of type II quasars. 

6. SUMMARY 

We have measured the 8 /xm luminosity func- 
tion of type I quasars as a function of redshift 
using 292 quasars and quasar candidates selected 
from Spitzer Space Telescope 24 /im imaging of 
the NDWFS Bootes field. Spectroscopic redshifts 
for 270 of these quasars were measured as part of 
the AGES survey of Bootes. We were able to de- 
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tect reddened type I quasars which were missed by 
UVX quasar surveys, and these comprise approx- 
imately a fifth of the type 1 quasar population. 
The type I quasar luminosity function at 8 /im, 
and its evolution with redshift, is very similar to 
quasar luminosity functions at optical and ultravi- 
olet wavelengths. Figure \W\ illustrates the evolu- 
tion of the space density of type I quasars. There 
is qualitative agreement between our best-fit lu- 
minosity function and the 2QZ and SDSS, despite 
our sample being selected at 24 iim rather than 
the optical. We directly detect the peak in quasar 
activity atz = 2.6±0.3, which has been difficult to 
directly measure with optical quasar surveys be- 
cause it corresponds to the redshift range where 
quasar optical colors are similar to those of stars. 

Our best-fit model indicates that the quasar 
space density declines by a factor of ~ 2 between 
z = 2.5 and z = 4.0, though this decline has 
low significance. If we fit our data with a model 
where the space density is constant after Zpeak 
then our luminosity function parameters change 
by < la. However, such a model disagrees with 
the results of optical surveys of 2; > 3 quasars 
(Schmidt, Schneider, & Gunn 1995; Fan et al. 
2004). Thus, unless the obscured quasar fraction 
increases at high redshift, our study and previ- 
ous optical quasar surveys indicate that the quasar 
space density declines at z > 3. 

While the evolution of the type I quasar space 
density has been measured over a broad range of 
wavelengths, the evolution of type II quasars re- 
mains uncertain. Figure [12] indicates that they 
could comprise 70% of z > 1 quasars brighter than 
1 mjy, though this fraction decreases with increas- 
ing flux and luminosity. Until there is a systematic 
campaign to measure the spectroscopic redshifts of 
type II quasars, the 8 /im luminosity function of 
all quasars will remain uncertain. 
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Table 1 

Compact and extended object number counts for the multiwavelength sample. 



24 jum Number counts 

flux (mjy) Total AGES Spectra R < 17.2 17.2 < R< 20.2 20.2 < ii < 21.7 21.7 < ii < 23 R> 23.^ 



Optically compact sources^ 
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23 


120 


67 


10 


2 


2 - 3 


54 


46 


5 


40 


9 








3 - 4 


21 


16 


5 


13 


3 








4 - 5 




6 


1 


5 


1 








5 - 6 


8 


3 


2 


4 


1 


1 





6 - 7 


6 


4 


1 


4 





1 





7 - 8 


4 


4 





4 











8 - 9 


3 


2 


1 


2 











9 - 10 























10 - 20 


9 


5 


5 


4 











20 - 30 


2 


1 


1 


1 











30 - 40 


1 





1 














40 - 50 























Optically extended sources 


1 - 50 


1965 


1181 


232 


1061 


296 


170 


206 


1 - 2 


1322 


735 


86 


700 


229 


138 


169 


2 - 3 


324 


224 


41 


201 


37 


17 


28 


3 - 4 


128 


87 


27 


74 


12 


11 


4 


4 - 5 


65 


43 


17 


33 


10 





5 


5 - 6 


33 


27 


10 


18 


3 


2 





6 - 7 


20 


13 


7 


9 


4 








7 - 8 


15 


14 


6 


9 











8 - 9 


9 


6 


4 


3 


1 


1 





9 - 10 


8 


6 


3 


5 











10 - 20 


24 


15 


15 


8 





1 





20 - 30 


12 


7 


12 














30 - 40 


3 


2 


3 














40 - 50 


2 


2 


1 


1 











^At faint magnitudes, 


where the SExtractor star-. 


galaxy classifier : 


is unreliable, we assume 


optical counterparts 


are extended. 


< 1 galaxies with AGES spectra that 


were classified as compact 


in the NDWFS imaging 


have been reclassified 


as extended objects. 



Table 2 

Luminosity Function Parameters 



Parameter 


All 


R - [24] > 4.4'' 


R - [24] < 4.4^ 


2QZb 


(A8(-26,0)[10-^Mpc"Vag-i] 
M-29, 2)[10-^Mpc-Vag-i] 


4.84±|:^| 


2.35 ± 0.28 


2.44 + 0.26 




4.46 ± 0.53 


2.16 ±0.26 


2.26 ±0.25 




a 


-2.75±0:1I 


-2.82 ±0.14 


-2.66 ±0.13 


-3.31 


ki 


+1.15 ±0.34 


+1.15 


+1.15 


+1.39 


k2 


-0.34 ±0.15 


-0.34 


-0.34 


-0.29 




±0.03 ±0.02 


+0.03 


+0.03 




^peak 


2.56 ±0.27 


2.56 


2.56 


2.4 


z range 


1 < z < 5 


1 < .2 < 5 


1< < 5 


0.4 < z < 2.1 


Spectroscopic sample size 


183 


81 


102 


15,830 



''Evolution parameters have been fixed for the red and blue subsamples. 

^2Q[L Zpeak value is an extrapolation of the evolution of the 0.4 < z <2.1 luminosity function. 
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Table 3 

Binned 1/Kiaa; z = 2 Quasar Luminosity Function 



8 yum absolute 


Space Density at z 


= 2 (Mpc-^mag-i)^ 


magnitude range 


1.0 < 2; < 5.0 Quasars 


1.5 < 2; < 2.5 Quasars 


-27.4 < Mg < -27.0 


4.3 ± 1.0 X 10"*^ 




-27.8 < Ms < -27.4 


2.4 ±0.5 X 10-*^ 




-28.2 < Ms < -27.8 


3.0 ±0.4 X 10-6 


2.4 ±0.6 X 10-*^ 


-28.6 < Ms < -28.2 


1.2 ±0.2 X 10"'^ 


1.1 ±0.2 X 10-6 


-29.0 < Ms < -28.6 


6.1 ±0.9 X 10"^ 


4.5 ± 1.0 X lO-'' 


-29.4 < Ms < -29.0 


2.6 ±0.5 X 10"^ 


2.7 ±0.7 X 10-^ 


-29.8 < Ms < -29.4 


2.4 ±0.5 X 10"^ 


2.0 ±0.6 X 10-^ 


-30.2 < Ms < -29.8 


9.2 ±2.4 X 10-8 


6.8 ±2.5 X 10-"^ 


-30.6 < Ms < -30.2 


4.7 ± 1.5 X 10-8 


7.4 ±2.8 X 10-8 


-31.0 < Ms < -30.6 


< 7.2 X 10^8 


< 2.8 X 10-"^ 


-31.4 < Ms < -31.0 


2.3 ±0.9 X 10-8 


4.9 ±2.1 X 10-8 



^Upper limits are 2>u. 
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